home *** CD-ROM | disk | FTP | other *** search
/ Monster Media 1996 #15 / Monster Media Number 15 (Monster Media)(July 1996).ISO / prog_c / cuj0696.zip / DWYER.ZIP / FREQ.TST / KSPROB.AWK < prev    next >
Text File  |  1995-09-30  |  649b  |  31 lines
  1. BEGIN\
  2. {
  3.     printf("Calculating Kolmogorov-Smirnov Probability K(N) <= X\n");
  4.     printf("Enter N: ");
  5.     Status = getline N;
  6.     if (Status <= 0)
  7.     {
  8.     exit;
  9.     }
  10.     N = N < 0 ? -int(N) : N == 0 ? N = 1 : int(N);
  11.     if (N < 0)
  12.     {
  13.     N = -N;
  14.     }
  15.     while (1)
  16.     {
  17.     printf("Enter Next X: ");
  18.     Status = getline X;
  19.     if (Status <= 0)
  20.     {
  21.         exit;
  22.     }
  23.     printf("For K(%d, %f):\n", N, X);
  24.     Darling = 1 - exp(-2*X^2)*(1 - (2/3)*X/sqrt(N));
  25.     Knuth   = 1 - exp(-2*(X + 1/(6*sqrt(N)))^2);
  26.  
  27.     printf("\tD. A. Darling's Approximation = %.14e\n", Darling);
  28.     printf("\tD. E. Knuth's   Approximation = %.14e\n", Knuth);
  29.     }
  30. }
  31.